Make Slow Fast - how to speed up interacting disordered matter 
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Anderson and dynamical localization have been experimentally observed with ultra-cold atomic 
matter. Feshbach resonances are used to efficiently control the strength of interactions between 
atoms. This allows to study the derealization effect of interactions for localized wave packets. The 
derealization processes are subdiffusive and slow, thereby limiting the quantitative experimental 
and numerical analysis. We propose an elegant solution of the problem by proper ramping the 
interaction strength in time. We demonstrate that subdiffusion is speeded up to normal diffusion 
for interacting disordered and kicked atomic systems. The door is open to test these theoretical 
results experimentally, and to attack similar computational quests in higher space dimensions 

PACS numbers: 05.45.-a, 71.55.Jv, 37.10.Jk 

Introduction 

The quantum wave nature of ultracold atoms in optical potentials [l[ , as demonstrated impressively through their 
macroscopic condensation 0, H| , is the key ingredient for the recent observation of Anderson localization with quantum 
atomic matter Hi]. Quasi -one-dimensional elongated traps are modulated randomly with speckle potentials [f| 0] , 
or simply quasiperiodically with interfering laser beams |8|, in order to observe the halt of spreading of an initially 
localized wave packet of 10 4 - 10 5 Rb and K atoms, and an exponentially localized atomic density distribution 
profile. The length scales are controlled by the localization length £ which is a function of the potential parameters, 
and possibly also the energies of packet atoms. This phenomenon of wave localization is inherently relying on the 
phase coherence of matter waves. It is closely related to the dynamical localization of the quantum kicked rotor in 
momentum space, which was successfully probed already in 1995 using ultra-cold Na atoms Q. Recent experiments 
with quasiperiodically kicked rotors with Cs atoms extend to two- and three-dimensional disorder potentials [loj . 
Interestingly systems of one- and two-dimensional optical waveguides have been also recently used to probe Anderson 
localization (III, H^. 

For some atomic species (KjCs, Na, Li) Feshbach resonances can be used to efficiently control the strength of 
interactions between atoms |L3l - [la | . This opens the possibility to study the fate of Anderson localization for interacting 
localized wave packets. Indeed the first experiment of this kind (l7| showed that interaction beats localization, but 
in a very slow way - the second moment mi of an atomic wave packet increases subdiffusively in time: ~ t a with 
a < 1. This process may stop in the long run once the atomic density n of the wavepacket reaches the inverse of 
the localization length £ (which touches the quantum world of many-body localization [l8j] ) . On shorter times (when 
typically more than 10 atoms occupy one local single particle state) the mean field approximation is a reasonable tool 
for the study of the subdiffusive process. 

The problem 

The mean-field approximation replaces the many body linear Schrodinger equation in a hugely dimensional Hilbert 
space with a nonlinear Schrodinger equation (NLS), e.g. the Gross-Pitaevsky equation. The effective interaction 
strength f3 is proportional to the scattering length a s . What matters in terms of quality here, is the fact that for almost 
any disorder (or quasiperiodic) potential realization the corresponding NLS will be nonintegrable. This seemingly 
unimportant mathematical property has a very profound impact - the dynamics of a wave packet becomes in general 
chaotic in time, characterized by positive Lyapunov coefficients and exponential divergence of nearby trajectories. As 
a consequence the coherence of phases of waves which constitute a given initial wave packet is lost, and with it also 
the whole effect of wave localization. First observed in 1993 by Shepelyansky for an NLS version of the quantum 
kicked rotor [li| , it was recently studied with great detail for random and quasiperiodic potentials [2(il - [24j . The main 
outcome for quasi-one-dimensional models with two-body atomic interactions is an asymptotic wavepacket spreading 
with a — 1/3 [25j |. With inverse time units equal to single particle kinetic energies the crossover from intermediate 
to universal asymptotic dynamics takes place at dimensionless time r ~ 10 6 . In the Florence setup [13] largest times 
reached are r ~ 10 4 , leaving the experiments in the intermediate case-to-case-dependent dynamics. Although the 
onset of subdiffusion is clearly observed, no reliable experimental data are currently at hand to measure the exponent 
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a, as follows from the data analaysis and the large statistical errors in [17] . While some experimental optimization 
and increase of the kinetic energy may add one order of magnitude in time, another one-two orders are needed and 
are probably currently out of range of accessibility. Notably similar problems of insufficient available time scales arise 
with computational studies when turning to higher dimensional analogs @, [2(| [2?J . While two-dimensional models 
appear to be at the edge of reasonable analysis, three-dimensional are clearly not. As follows from the above we can 
diagnose the problem of lacking time scales for a safe observation and study of subdiffusive interacting atomic cloud 
dynamics in disordered media. 

The solution 

Instead of trying to substantially increase available time scales, we propose here to speed up the subdiffusive process 
itself. This is done by a temporal ramping of the two-body interaction strength, which can be varied e.g. for K atoms 
by three orders of magnitude close to the Feshbach resonance [l3[ . Why should that help? The momentary diffusion 
rate D of a spreading packet in one spatial dimension is proportional to the fourth power of the product of interaction 
strength (3 and particle density n: D ~ (fin) 4 (28J. In the course of cloud spreading the density n decreases, and 
therefore also D. This is the reason for the predicted subdiffusion process, which is substantially slower than normal 
diffusion. We propose here to compensate the decrease of the density n with an increase in the interaction strength 
(3. Depending on the concrete time dependence f3(r) we expect different faster subdiffusion processes, and possibly 
even normal diffusion. The condition for that outcome to be realized is, that the internal chaos time scales (basically 
the inverse Lyapunov coefficients) will be still short enough so that the atomic cloud can first get chaotic, and then 
spread. With that achieved, the cloud spreading will be faster, and we can expect that the available experimental 
time will suffice for the precise observation and analysis of the process. 

Let us get into numbers for one spatial dimension. The second moment is mi ~ 1 /n 2 and the momentary diffusion 
constant D ~ (/?n) 4 . For a constant (3 the solution of mi — Dt yields mi ~ 1/n 2 ~ r 1 / 3 , and therefore n ~ r -1 / 6 . 
Thus we choose now a time dependence /3 ~ t v . Then the resulting spreading is characterized by 

m2 „ T (.l+*»)/3 , d = 1 . (1) 

For v = 1/2 we already obtain normal diffusion mi ~ r. 

Similar for two spatial dimensions, where m 2 ~ 1/n, for a constant /3 the cloud spreading is even slower with 
mi ~ r 1 / 5 . With a time dependent ramping /3 ~ t v the resulting speedup is 

m 2 ~ r (l+ 4 -)/ 5 , d = 2 . (2) 

For v =1 we again obtain normal diffusion. 

Once ramping is too fast, we expect to see several different scenaria. Either fragmenting atomic clouds appear 
since some parts of the cloud get self-trapped [29|, [3(| and some other parts do not. If self-trapping is avoided, we 
may also see ramping-induced diffusion: while the internal cloud dynamics does not suffice to decohere phases, initial 
fluctuations in the density distribution can lead to considerably different temporal energy renormalizations in different 
cloud spots, and therefore to an effective dephasing similar to a random noise process in real time and space. 

Results in one dimension 

Here we study the spreading of atomic clouds in one-dimensional disorder potentials and in a quantum kicked rotor 
with interacting atoms. The first model is described with the discrete NLS (DNLS) 

.din _ ^ + p ^ |^|2 ^ _ ^ _ ^ ^ 
or 

in which the on-site energy e; is chosen uniformly from a [-W/2, W/2] random distribution. The nonlinear quantum 
kicked rotor (NQKR) is studied within the diagonal interaction approximation introduced by Shepelyansky in [19] : 

Mr + 1) = ^(- i )'-" l J i _ m (fc)V^(r)e- I 5 m2+I ' 3 ^)l^l 2 , (4) 
m 

where 4>i(t) are the Fourier coefficients of the corresponding time-dependent many body wave function. Ji_ m (k) is a 
Bessel function of the first kind, whose argument k is the kick strength, and f is a parameter which relates the period 
of applied kicks T (set to T = 1) to the natural frequency of rotor, defined as cj = H/2M (M-mass of atoms). In both 
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models is the interaction strength ramped in time r - the dimensionless time for the DNLS and the number of kicks 
for the NQKR model: 

(00 , T < T 

«"-{*(*)" ,r>r a < 5 > 

In both models we consider a wave packet which is initially concentrated on a sin gle s ite for purely technical reasons, 
without any loss of generality (extended initial clouds are perfectly usable as well [3l|). After some first time scale To 
the packet spreads approximately over one localization length £. The total norm of the packet is set to one without 
any loss of generality and is proportional to the total number of atoms in a cloud (similarly, we could also choose 
any larger norm and rescale j3 accordingly). To characterize the spreading of the cloud we compute the density 
n i = I'M 7 ") I ■ the participation number P = n f (the- number of strongly excited sites), and the second moment 

m-2 = J^r(Z — l) 2 ni (the squared distance between the wave packet tails), where I = Y]j Irii is the first moment. In 
the NQKR the average energy of the atomic cloud is proportional to the corresponding second moment, E = \mi. 
Another remarkable difference between both systems is that for large values of self-trapping can occur for atomic 
clouds in disordered spatial systems which may lead to soliton formation 32]. For the kicked rotor case this is 
impossible since the interaction strength ft action in (U) is cyclic reflecting the circumstance of periodic kick action in 
momentum space [33| . 

Equation (j3|) was time evolved using a SABA-class symplectic integration scheme and equation (j4|) as an iteration 
map. The parameters were fixed to (3q = 1 and W = 4 for the DNLS and /?o = 0.4, f = 1 and k — 3 for the 
one-dimensional NQKR. Different realizations for the DNLS were produced by choosing different unique random 
sequences in the interval [-W/2, W/2], while for the NQKR they were realized by exciting different initial states. 
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FIG. 1: Evolution of the averaged norm density < ni(r) > in the case without (y = 0) and with ramping {y = 0.3) in log scale 
for the DNLS model. 

The spreading of wave packets in the DNLS model, without and with ramping of the nonlinearity are shown in Fig. 
(fTj). Clearly packets spread faster when the nonlinearity is ramped in time. To quantify the spreading exponent, we 
averaged the logs (base 10) of P and over 1000 different realizations and smoothened additionally with locally 
weighted regression [34j]. The (time-dependent) spreading exponents are obtained through central finite difference 
method [Hi, a = d ■ The results for the DNLS and NQKR model are shown in Fig. ©. The exponents 
of subdiffusive spreading reach the theoretically predicted values. Moreover, for faster ramping of nonlinearity, the 
asymptotic state with constant exponent is reached faster. Monitoring of the participation number P for the DNLS 
indicates that self-trapping starts to occur already for v — 0.4. Results for the NQKR model, in which the self-trapping 
is avoided, confirm the reaching of a normal diffusion process for v = 0.5. Remarkably, the absence of self-trapping 
for the NQKR results in superdiffusion on intermediate times for v > 0.5. Finally, the exponent relaxes back to the 
normal diffusion value indicating the realization of ramping-induced diffusion. This case is illustrated for v — 1.5 in 
Fig. ©. 
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FIG. 2: Left column: the second moments (upper) and their power-law exponents a (lower) for the DNLS model for v = 
(red), v = 0.1 (green), v — 0.2 (blue), v — 0.3 (magenta), v — 0.4 (cyan), and v = 0.5 (black). Right column: the second 
moments (upper) and their power-law exponents a (lower) for the NQKR model for v = (red), v — 0.17 (green), ^ = 0.25 
(blue), i/ = 0.33 (magenta), v — 0.5 (cyan), and v = 1.5 (black). Dashed colored lines correspond to expected values for 
exponents in both cases. 



Results in two dimensions 

To speed up subdiffusive processes in higher dimensions we considered the two-dimensional NQKR model based 
on the map introduced in [361 ] . with an additional phase term which takes into account interactions in the diagonal 
approximation: 

VW 2 (t + 1) = ^(-*) Sl+s V Sl (fc/2)J S2 (fc/2)^ 1 _ Sl _ S2 , /2+Sl _ S2 (r) 

si,s 2 

e -i|((ii-si-s 2 ) 2 + (i2+si-s 2 ) 2 )+i0(T)|'0 il _ sl _ 32>i2 + 31 _ 32 (T)| 2 ^ fQ\ 

The notation is the same as in the one-dimensional case, except that now we have two indices, for each possible 
direction. Note that according to the relation for Bessel functions J- n {x) — (—l) n J n (x), the wave packet ifii lt i 2 , defined 
by expression ([6]), exhibits symmetry with respect to and I2 direction. The density is defined as fi^ j a = |V , ii,/ 2 ( T )l > 
the participation number as P = 1/53/ i 2 n h / 2 > an< ^ the secon d moment is 777,2 = J2h i 2 [(h — + — h) 2 }ni lt i 2 , 
where l\ — i 2 l\ni lt i 2 and I2 = 53; i 2 l2 n h,h- Equation ([6]) was time evolved as an iteration map for the fixed 
parameters (3q = 0.4, f = 1 and k = 2. 

In Fig. ((3]) we compare the wave packet evolution for v — 1 (normal diffusion) and v = at three different moments 
of time. We clearly observe the symmetry of wave packet and a much more violent spreading in the presence of 
ramping. The spreading exponents are computed similar to the one-dimensional case. We find very good agreement 
with the theoretical prediction (Fig. (g])). Again the asymptotic spreading state is reached faster for stronger ramping. 

Conclusion 

We have investigated the speeding up of the subdiffusive spreading in interacting disordered and kicked atomic 
systems by a proper ramping of the interaction strength in time. We confirm that ramping the interaction strength 
leads to faster subdiffusion. For fast enough ramping we even reach normal diffusion of atomic clouds. Self-trapping 
effects in disordered systems are limiting further speed up of the wave packet spreading. Most importantly the concept 
works equally well in one-dimensional and two-dimensional systems. 

Our results on how to speed up slow subdiffusive processes in interacting disordered matter will be useful for 
quantitative experimental and computational studies of the impact of interactions on disorder induced matter wave 
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FIG. 3: The norm denisty ni lt i 2 (r) in the case without (y = 0) (upper row) and with ramping (y — 1.0) (lower row) after 
t — 10 2 , r = 10 3 , and r = 10 5 kicks in log scale for the two-dimensional NQKR model. 
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FIG. 4: The second moments (upper) and their power-law exponents a (lower) for the two-dimensional NQKR model for v = 
(red), v — 0.5 (green), v — 0.75 (blue), and v = 1 (magenta). Dashed colored lines correspond to expected values for exponents. 



localization. This is particularly true for experimental realizations with ultra-cold atoms, where the scattering length 
and thus the interactions strength can be tuned via Feshbach resonances by magnetic field variations. On the other 
hand, our results are also very useful for computational studies of spreading regimes in higher dimensional systems, 
where even modern computers reach their limits before reaching the subdiffusive asymptotics. 
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